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NOMENCLATURE 


a design variables, eqs. (1) 

b design variables, eqs. (1) 

c airfoil chord, m (ft) 

c, section drag coefficient 

d 

C£ section lift coefficient 

c section pitching-moment coefficient referenced to quarter chord 
m 

Pi - Poo 

C pressure coefficient 

P ‘loo 

C * pressure coefficient for M = 1 
P * 

f shape function, eqs. (1) 

g shape function, eqs. (1) 

h tunnel height, m(ft) 

M free-stream Mach number 

m shape function exponent, eqs. (2) 

n shape function exponent, eqs. (2) 

p static pressure, N/m^ (Ib/ft^) 

q dynamic pressure, N/m^ (Ib/ft^) 

R Reynolds number based on free-stream conditions and airfoil chord 
r shape fimction exponent, eqs. (2) 

s shape f motion exponent, eqs. (2) 

X airfoil abscissa, m (ft) 

y airfoil ordinate, m (ft) 

Subscrip ts : 
max maximum 

local ' 


iii 



free stream 


00 


6 wall correction factor 

o angle of attack 
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SUMMARY 


A 10-percent thick helicopter rotor section designed by numerical optimi- 
zation has been tested at Mach numbers of 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.75, 
0.77, 0.79, 0.8, 0.81, 0.82, 0.83, and 0.84 with Reynolds numbers ranging from 
1.9x10® at Mach 0.2 to 4x10® at Mach numbers above 0.5. During the design 
process, emphasis was placed upon achieving acceptable maximum lift coeffi- 
cients at Mach numbers below 0.6, a drag divergence Mach number of at least 
0.8 at lift coefficients near 0 and low pitching-moment coefficients. 

Test results have shown that the maximum lift coefficient decreased from 
approximately 1.33 at Mach numbers of 0.2, 0.3, and 0.4 to below 1 at Mach 0.7. 
A drag divergence Mach number of 0.82 at lift coefficients near 0 was achieved 
with some "drag creep" at Mach numbers above 0.6. Pitching-moment coefficients 
at zero lift decreased from approximately -0.01 at Mach 0.2 to -0.03 at 
Mach 0.8. 


INTRODUCTION 


An automated airfoil section design program utilizing a conjugate gradi- 
ent numerical optimization algorithm (ref. 1) coupled with an aerodynamic 
analysis program (ref. 2) has been used successfully in the past to design 
advanced airfoil sections at low to transonic speeds (refs. 3-5). A major 
advantage of design by numerical optimization is that the designer may consi- 
der multiple design points simultaneously during the design process. For 
example, "drag creep" was eliminated during the design of an advanced super- 
critical airfoil section for Mach 0.78 by imposing a drag constraint at Mach 
0.76 (ref. 3). The constraint prevented the optimization code from seeking a 
profile which was highly design-point optimized without regard to off-design 
performance. Such direct consideration of two design conditions would be 
difficult, if not impossible, with inverse design methods. This feature of 
design by numerical optimization makes it useful in the design of rotor sec- 
tions since several design conditions must be considered simultaneously dur--'^" 
ing the design process. Some of these design conditions are high maximum-lift 
coefficients and good stall characteristics from Mach 0.3 to 0.5, a high 
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/ rag ratio at Mach 0.6 and Cj^ ~ 0.6, a drag divergence Mach number of at 
east 0.8, no drag creep , and low pitching moments over most of the Mach 
number range. To achieve all these conditions the optimization code must 
monitor and constrain certain aerodynamic parameters at five or six different 
combinations of Mach number and angle of attack. Such a design problem 
inevitably means compromise and would normally be carried out by trial and 
error with the designer making all decisions and providing all computer mani- 
pulations. With numerical optimization, the computer makes many of the design 
decisions thus reducing the time and manpower required during the design pro- 
cess. While there are no fundamental limits on the number of design points 
that can be considered, there are cost limitations. For example, a typical 

problem would use approximately 8 hr of CPU time on a CDC 
7600 computer. -During the study reported here, three design conditions were 
considered. Two conditions addressed the high-lift performance of the retreat- 
ing blade while the third considered advancing blade performance. "Drag-creep" 
and hover performance were not considered in order to reduce computing time 
and cost. As will be shown later, omitting "drag creep" from consideration 
compromised airfoil performance above Mach 0.7. 


DESIGN METHOD 


^ Only a brief description of numerical optimization will be given here 
since a complete discussion of the technique can be found in reference 6. 

Airfoil design by numerical optimization starts with a baseline profile 
provided by the designer. The baseline profile can be somewhat arbitrary but 
is usually chosen from a family of airfoils which exhibit some of the charac- 
teristics desired by the designer. Such a choice reduces the time and cost 
required to finalize the design. 

The airfoil geometry is represented in the design program by the follow- 
ing equations 
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where yus^gggj^^j^g, ^^^baseline ordinates of the upper and lower 

surfaces of the baseline airfoil, a^ and b^ are the participation coefficients 
(design variables) and f^ and g^j^ are the shape functions which are added 
linearly to the baseline profile by the optimization program to achieve the 
desired design Improvement. The contribution of each function, f^ or g^, is 
determined by the value of the participation coefficients a^ and bi. The 
shape functions, f^ and gi, can be any smooth function defined over the air- 
foil chord, but experience has shown that the following functions give accept- 
able results from most design problems: 
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or = sin[ir(x/c) ] 


or 


f „ (x/c)®[l - (x/c)] 

^1 Si ToJTJ) 


C2) 


where n, m, s, and r are exponents which determine the chordwise locations 
and/or magnitude of each shape function. Such functions are well suited to 
airfoil design because the effect of each function can be concentrated over a 
limited region of the airfoil. 

The optimization design process begins by setting all participation coef- 
ficients to zero so the first computation gives the aerodynamic characteristics 
of the baseline profile. After computation of the characteristics of the base- 
line profile the optimization program changes the participation coefficients, 
one by one, from the initial value of 0 to a small number e << 1. The aero- 
dynamic analysis of each perturbed airfoil was obtained and used to calculate 
the gradient of the aerodynamic coefficients with respect to each design vari- 
able. The importance of each shape function in achieving the design goals can 
be ascertained at any step during the optimization process by observing the 
value of the participation coefficient associated with each function at that 
step. 


The optimization design program used during this study permits the 
designer to designate one or more of the aerodynamic quantities to be mini- 
mized or maximized as the objective function. All other aerodynamic and geo- 
metric quantities pertinent to the design are treated as constraints. The 
designer must decide which parameters will be treated as constraints and which 
will be the objective function. If the choice is not obvious it will be 
necessary to use trial and error to determine which combination of objective 
function and constraints gives the best overall result. 

During the present study, the shock-drag coefficient at M = 0.82, a = 0® 
was chosen to be the objective function. Constraints were imposed on the 
shock-drag coefficients (cj < 0.001) at M = 0.4, a = 12®, and M = 0.5, a = 10® 
to delay retreating blade shock-stall. The shock-drag objective function 
forced the optimization program to seek a shock-free airfoil at M = 0.82, 
a = 0®. As the optimization process proceeded it became clear that the pro- 
file required for good advancing blade performance would not necessarily pro- 
duce good retreating blade performance. Hence, the designer had to decide if 
the drag constraints imposed at M = 0.4 and 0.5 were too restrictive and 
whether or not to require a shock-free profile for M = 0.82. In other words, 
the design required a compromise and it was decided that the design process 
should not be completely automated for this first application of numerical 
optimization to the design of a rotor section. Hence, the designer remained 
in the "loop" for the entire optimization process monitoring design progress, 
adjusting shape functions, altering constraints, and changing design conditions 
in an effort to find the most suitable compromise. 
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The aerodynamic analysis code contained in the design program does not 
consider the effect of viscosity. Hence, several intermediate airfoils were 
selected during the design process for boundary layer and stall analysis. If 
boundary-layer growth appeared excessive or premature separation was indicated, 
adjustments were made to the shape functions and/or constraints to guide the 
optimization process to a more efficient airfoil. 

The airfoil selected as the initial profile to start the optimization 
process for this design problem was a Wortman airfoil designed for helicopter 
applications. This profile was chosen since it already exhibited some of the 
desired characteristics of a rotor section and hence would require less com- 
puter time to achieve the final design. The profile resulting from this 
demonstration of design by numerical optimization is designated the A-1 air- 
foil since it is the first rotor section designed by numerical optimization at 
Ames Research Center. 


APPARATUS AND TEST PROCEDURE 


Model 

An airfoil model of the A-1 airfoil (shown in fig. 1) was machined from 
a stainless-steel billet (the profile coordinates are given in table 1). The 
model had a chord of 15.24 cm (6 in.) and a span of 60.96 cm (24 in.). The 
model was equipped with 24 upper-surface orifices and 23 lower-surface orifices 
drilled normal to the surface to determine the pressure distribution on the 
model surfaces. 


Wind Tunnel 

The test was conducted in the Ames 2- by 2-Foot Transonic Wind Tunnel, a 
variable-speed, continuous-flow, ventilated-wall , variable-pressure facility . 
The tunnel can be used for two-dimensional testing by replacing the ventilated 
side walls with solid walls where model-supporting thick glass windows are 
mounted. The windows are rotated by a motorized drive system to change the 
angle of attack. A 82-tube drag rake located 1.75 chords behind the model 
trailing edge was used to survey the model wake. Airfoil models are mounted 
spanning the horizontal dimension of the tunnel test section so that the 
center bf rotation is near the 25-percent chord station on the model. The 
gaps between the side walls and the end of the model were sealed. 


Instrumentation 

Measurements of the model surface pressures and the wake-rake pressures 
were'* made by an automatic pressure-scanning system that utilizes precision 
pressure transducers. Basic tunnel pressures were measured with precision 
mercury manometers. Angle of attack was measured with a potentiometer operated 
by the drive gear for the rotating side windows. Data were obtained by a 
high-speed, data-acquisition system and recorded on paper tape. 
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Tests 


The section aerodynamic coefficients of the A-1 airfoil model were 
obtained at the following nominal test conditions: 


M Rx10“6 


0.20 

1.9 

.30 

2.5 

.40 

3.5 

.50 

3.5 

.60 

4.0 

.70 

4.0 

.75 

4.0 

.77 

4.0 

.79 

4.0 

.80 

4.0 

.81 

4.0 

.82 

4.0 

.83 

4.0 

.84 

4.0 


The angles of attack ranged from approximately - 2 ° to stall. At the higher 
Mach numbers it was not possible to obtain data for a fully-stalled condition 
because of severe vibration of the wake rake due to the unsteady loads 
induced by total flow separation over the upper surface of the model. Data 
were obtained for all test conditions with free transition only. 

Pressure coefficients were determined from surface pressure measurements. 
Section normal-force coefficients, chord-force coefficients, and pitching- 
moment coefficients were obtained from an integration of the pressure coeffi- 
cients. The pitching-moment coefficients were referenced to the quarter-chord 
point- Section profile drag was calculated from the wake-rake total and 
static-pressure measurements. 

The model angle of attack was corrected for the presence of the tunnel 
walls by the following equation: 


Aa = 6(c/h)c^ 

where Aa, 6, c/h, and C£ are the angle-of-attack correction, correction 
factor, model chord/tunnel height ratio, and section lift coefficient, respec- 
tively. The angle-of-attack correction factor 6 is a function of Mach 
number. The following values were used and the corresponding Aa was added 
algebraically to the model geometric angle of attack expressed in degrees: 
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M 

ol 

0.2 

5.4 

.3 

8.5 

.4 

10.6 

.5 

11.9 

.6 

12.5 

.7 

12.8 

.75 

12.9 

.77 

12.9 

.79 

12.9 

.80 

13.0 

.81 

13.0 

.82 

13.0 

.83 

13.1 

.84 

13.1 


RESULT AND DISCUSSION 


The basic aerodynamic force coefficients for the A-1 airfoil are shown 
in figures 2(a) through 2(n). The lift curves indicate that the stall charac- 
teristics are fairly abrupt at Mach numbers 0.2, 0.3, and 0.4 but become much 
more gradual at Mach numbers greater than 0.4 (compare fig. 2(a) and 2(e)). 

The stall characteristics were difficult to determine at the higher Mach num- 
bers due to severe vibration of the wake rake under conditions of massive flow 
separation. The lift curves exhibit some nonlinearity at the lower Mach num- 
bers which may be due to variations in transition location. 

The drag polars at Mach numbers below 0.7 show a minimum drag point at 
lift coefficients above 0.3 due to an increase in laminar flow over the lower 
surface which is caused by the pressure gradient changing from adverse to 
favorable near the leading edge (compare the pressure distributions in figs. 
3(a) and 3(d)). 

The pitching— moment curves exhibit some nonlinearity with increasing lift 
coefficient which may be attributed to variation in the transition location. 
The pitching-moment coefficient at zero lift decreased from approximately 
-0.01 at Mach 0.2 to -0.03 at Mach 0.8. Pitching moments of this magnitude 
are somewhat greater than desired and are due primarily to the inability of 
the aerodynamic theory contained within the design program to predict the 
pressure distribution accurately near the trailing edge. 

Pressure distributions for all test conditions are shown in figures 3 
through 16. Of particular interest is the pressure "peak" appearing near the 
leading edge on the lower surface at low lift coefficients for all test Mach 
numbers. This pressure "peak" is due to the flow accelerating too rapidly 
around the leading edge of the lower surface and is a result of too much 


iThe correction factors, 6, were determined during a tunnel calibration 
conducted by L. S. Stirers, Jr. 
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forward curvature. The pressure "peak" could have been minimized during design 
by imposing a constraint on shock drag at a Mach number slightly less than 0.8 
to avoid "drag creep." Such a constraint would have forced the optimization 
code to seek a curvature distribution which would accelerate the flow more 
gradually along the lower surface, thereby delaying shock development until 
after "drag rise." The pressure distributions along with the drag curves 
clearly show that the lower surface pressure "peak" was the main contributor 
to the "drag creep" of the A-1 airfoil. Further work is needed to assess the 
effect of changing the leading-edge contour on the overall performance of the 
profile. 

The upper surface pressure distribution exhibits a slight pressure "hiimp" 
over the forward region of the profile at low lift coefficients which gives 
rise to a weak shock at transonic conditions. The moderate adverse pressure 
gradient following this "hump" would cause transition to occur near the lead- 
ing edge at low Mach numbers but would not cause separation for typical flight 
Reynolds numbers as indicated by the good trailing-edge pressure recovery. 

Summary curves of drag coefficient vs Mach number for lift coefficients 
from 0 to 0.5 are shown in figure 17. The "drag creep" due to the pressure 
"peak" near the leading edge on the lower surface discussed earlier is clearly 
evident above Mach 0.6 at c^ = 0. As the lift coefficient increases, the 
pressure "peak" on the lower surface decreases and the "drag creep" is replaced 
by an early "drag rise" due to increasing shock strength on the upper surface. 

Summary curves of pitching-moment coefficient vs Mach number for lift 
coefficients from 0 to 0.5 are shown in figure 18. The desired range of 
-0.01 < Cjjj < 0.01 was achieved only at Mach 0.2 for all lift coefficients but 
was attained at Mach numbers to 0.7 at' lift coefficients of 0.4 and 0.5. The ' 
difficulty in achieving low pitching moments was due primarily to two factors: 
first, the aerodynamic theory contained in the design program did not predict 
the pressures near the trailing edge well enough, and second, the magnitude 
of the pressure "peak" near the leading edge on the lower surface was under- 
predicted by the theory for some flight conditions causing the importance of 
the "peak" to be underestimated. If the "peak" were eliminated the suction 
pressures over the forward lower surface would be minimized thereby reducing ' 
the nosedown pitching moment. 

A summary plot of maximum lift coefficient vs Mach number is shown in 
figure 19. The maximum lift coefficient remains nearly constant from Mach 0.2 
to Mach 0.4 at a value of approximately 1.33 and falls rapidly thereafter. 

The constancy of maximum lift coefficient at lower Mach nianbers is due in part 
to increasing Reynolds number with increasing Mach number, but it also reflects 
the favorable pressure distributions generated by the design method at 
M = 0.4^ 

Some typical experiment-theory comparisons are shown in figure 20. In 
general, the theory predicts the basic character of the flow fairly accurately. 
However, there are discrepancies between experiment and theory in some areas. ' 
First, the trailing-edge pressure recovery is consistently overestimated "by' 
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theory, even for attached flow. Second, the shock location shown by the exper- 
imental data is somewhat further aft than indicated by theory (figs. 20(g) and 
20(h)). Finally the theory tends to give poor drag prediction. A rather 
surprising feature of the experiment-theory comparison is the good agreement 
for the low Mach number, high lift condition (e.g., fig. 20(b)). This agree- 
ment was not anticipated because the Nash— McDonald, integral, boundary— layer 
method used to correct the potential aerodynamic theory was developed for 
flattop transonic pressure distributions with concave pressure recoveries. 

The fact that the theory reproduces the main features of the flow reasonably 
well makes a redesign effort of the A— 1 airfoil possible. A complete descrip- 
tion of the theoretical aerodynamic analysis program used for the experiment- 
theory comparison can be found in reference 7. 


CONCLUDING REMARKS 


An advanced helicopter rotor airfoil section has been designed by numer- 
ical optimization. The section is intended for use on rotor blades of con- 
stant profile and therefore was designed to exhibit acceptable aerodynamic 
characteristics over a wide range of Mach number, angle— of— attack conditions. 
Such a section would be expected to be less than optimum at each condition. 

-.Wind-tunnel test results have shown that a drag divergence Mach number 
of 0.82 at zero lift was achieved with a maximum lift coefficient of approxi- 
mately 1.33 at Mach numbers of 0.2, 0.3, and 0.4. Test Reynolds numbers varied 
from 1.9x10® at Mach 0.2 to 4x10® for Mach numbers greater than or equal to 
0.6. A modest amount of "drag creep" occurred at low lift coefficients at 
Mach numbers above 0.6 due to rapid acceleration of the flow to supersonic 
speed along the lower surface. The leading-edge bluntness which caused the 
acceleration was developed to promote acceptable stall characteristics and 
C£jnax Mach numbers. More work is needed to find a better 

compromise profile which will exhibit good "drag-rise" characteristics, accept- 
able "drag creep." 

In spite of the deficiencies found with the A-1 airfoil, numerical opti- 
mization has been found to be a good technique for rotor section design; the 
main advantage being the ability to consider multiple design conditions simul- 
taneously during the design process. 

The test results validated the basic approach for the multiple design 
point problem of rotor blades with conflicting requirements. 
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TABLE 1 .- COORDINATES FOR THE 
A-1 AIRFOIL 


0.00000 

0.00000 

0.00000 

.00020 

.00238 

-.00233 

.00050 

.00377 

-.00338 

.00100 

.00541 

-.00472 

.00200 

.00766 

-.00651 

.00350 

.01013 

-.00844 

.00500 

.01214 

-.00995 

.00650 

.01388 

-.01120 

.00800 

.01543 

-.01227 

.01000 

.01732 

-.01350 

.01250 

.01945 

-.01482 

.01600 

.02214 

-.01634 

.02000 

.02490 

-.01777 

.02500 

.02801 

-.01922 

.03500 

.03335 

-.02137 

.05000 

.03991 

-.02365 

.06500 

.04523 

-.02549 

.08000 

.04961 

-.02710 

.10000 

.05421 

-.02902 

.12500 

.05829 

-.03104 

.15000 

.06098 

-.03277 

.20000 

.06344 

-.03551 

.25000 

.06431 

-.03727 

.30000 

.06446 

-.03828 

.35000 

.06409 

-.03866 

.40000 

.06316 

-.03848 

.45000 

.06154 

-.03782 

.50000 

.05924 

-.03665 

.55000 

.05623 

-.03501 

. 60000 

.05249 

-.03297 

.65000 

.04792 

-.03056 

.70000 

.04246 

-.02785 

.75000 

.03600 

-.02486 

.80000 

.02860 

-.02153 

.85000 

.02064 

-.01786 

.90000 

.01260 

-.01374 

.92500 

.00899 

-.01144 

.95000 

.00598 

-.00888 

.97500 

.00392 

-.00603 

.99000 

.00322 

-.00421 

1.00000 

.00299 

-.00300 
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(a) Concluded. 
Figure 2 .- Continued. 
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(b) - 0.3, - 2.5i<10‘. 

Figure 2.- Continued. 
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(b) Concluded. 
Figure 2.- Continued. 
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(c) Concluded. 
Figure 2.- Continued. 
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Figure 2.- Continued. 
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(e) Concluded. 
Figure 2.- Continued. 
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Figure 2.- Continued. 
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(h) M 


set 


= 0.77, Re 


set 


4x10®. 
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Figure 2.- Continued. 
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Figure 3.- Pressure distributions for A-1 airfoil, = 0.2. 
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Figure 3.- Continued. 
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Figure 3.- Continued. 
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Figure 3.- Continued. 
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Figure 4.- Pressure distributions for A-1 airfoil, 
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Figure 8.- Pressure distributions for A-1 airfoil, M =0.7. 
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Figure 9.- Pressure distributions for A-1 airfoil, *» 
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Figure 10.- Pressure distributions forA-1 airfoil, M 
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Figure 12.- Pressure distributions for A-1 airfoil, *= 0.8. 
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Figure 13.- Pressure distributions for A-1 airfoil, = 0.81. 
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Figure 14.- Pressure distribution for A-1 airfoil, M = 0.82. 
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Figure 16.- Pressure distributions for A-1 airfoil, M = 0.84. 
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Figure 19.- Variation of section maximum- lift coefficient with Mach number 

for the A-1 airfoil, transition free. 
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